Spin gap in low-dimensional Mott insulators with orbital degeneracy 
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We consider the exchange Hamiltonian Hst = — J X/<rr'> (2S r ■ S r ; — |)(2T r ■ T r i — 5) , describing 
two isotropic spin-1/2 Heisenberg antiferromagnets coupled by a quartic term on equivalent bonds. 
The model is relevant for systems with orbital degeneracy and strong electron-vibron coupling in 
the large Hubbard repulsion limit. To investigate the ground state properties we use a Green's 
Function Monte Carlo, calculating energy gaps and correlation functions, the latter through the 
forward walking technique. In one dimension we find that the ground state is a "crystal" of valence 
bond dimers. In two dimensions, the spin gap appears to remain finite in the thermodynamic limit, 
and, consistently, the staggered magnetization - signal of Neel long range order - seems to vanish. 
From the analysis of dimer-dimer correlation functions, however, we find no sign of a valence bond 
crystal. A spin liquid appears as a plausible scenario compatible with our findings. 



The object of this study is a particular spin-exchange 
Hamiltonian in which two spin-1/2 variables - a stan- 
dard spin S, and a pseudo-spin T representing an orbital 
degree of freedom - are coupled together in the following 
way: 
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the lattice summation running over the nearest neighbor 
sites of a lattice. The symmetry group displayed by Eq. 
(Q) includes SU(A)E This Hamiltonian was shown to de- 
scribe the low-energy physics of an insulating crystal with 
one-electron per site in a two-fold degenerate orbital, in 
the limit of large on-site repulsion (Mott insulator) and in 
presence of Jahn- Teller effect.El The derivation of-Ea. M), 
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henceforth referred to as ST-model, is standardErta The 
crucial physical condition to be verified is that, among 
the possible two-particle states obtained upon virtual 
hopping, the inter-orbital singlet should turn out to be 
the lowest in energy, which is indeed the case when a 
strong dynamical JT effect is at play. A different (per- 
haps more common) physical situation is that the lowest 
two-particle intermediate state is a triplet (as a result of 
a Hund's rule), in which case the exchange Hamiltonian 
has a differentpform and spin ferromagnetism is the natu- 
ral outcome In the general case, the exchange Hamil- 
tonian contains pseudo-spin anisotropic terms, and its 
phase diagram-has been addressed with different tech- 
niques in Refs.0i3 n 
Ceo charge transfer compounds, based on Cq ions,u 
or molecular compounds with similar characteristics, are 
potential candidates to the realization of a scenario in 
which the Hamiltonian in Eq. (Q) is possibly relevant. 
Cgo has a triply degenerate t\ u molecular orbital, coupled 
to several intramolecular vibrations, resulting in an im- 
portant dynamical JT effect .13 The electron-phonon cou- 
pling leads to a substantial pairing energy m Q.leV which 
is, however, overwon by a substantially larger Hubbard 
PrjI - 1.5eV,EI resulting in a Mott insulating behavior. 
The lattices tend to have anisotropic lattice constants, 



with a pronounced quasi-one-dimensional character B As 
for the virtual intermediate states of Cqq , singlets and 
triplet are close in energy, at least in solutions .El Ifr)^ 
triplet prevails, as is possibly the case in TDAE-C6oP : t2l 
there can be spin ferromagnetism. We believe, however, 
there is room in some future system for the alternative 
possibility of a spin gap, and the formation of the spin- 
orbital VB state described here. 

Simple variational arguments are quite instructive in 
elucidating, as a function of the dimensionality of the 
lattice, the physical role played by the quadratic and 
quartic terms in which the ST-model can be rewritten, 

HsT = J J2<rr'>{Sr'Sr' +T r -T r t — 4(S r "S r ') (T r -T r /)J. 

The quadratic terms are standard antiferromagnetic cou- 
plings, which, alone, would drive the system to a Neel 
ordered state for both spin and orbital variables. The 
quartic term, on the contrary, gains from large values 
of (S r • S r /) and (T r • T r /): the optimal situation is a 
singlet valence bond (VB) state for both S and T vari- 
ables, where (S r ■ S r >) = (T r ■ T r >) = -3/4. Let us 
compare the energy of the following two states: (i) The 
product of two independent Heisenberg ground states, 
|*2-h) = I^h )®!*^), and (ii) the product oftwoVB 

states, |* 2 -vb) = \^vb) ® |*vb) where each VB state 
has the simple form of a product of dimers on adjacent 
sites. In D = 1 we take |* V b) = (12) (34) ■■■ (L-1L), 
where (ij) denotes a singlet between sites i and j. In 
D = 2 we consider any state with nearest-neighbor 
pairs coupled into singlets. The energy contribution 
of the quartic term is substantially lower for the VB 
state, — 4J(l/2)(— 3/4) 2 , than for the Heisenberg state 

for D < 2, -4JD(e n /DJ) 2 , where e H /J = £>(^ S) |S r • 



S r ,|* 



-0.4431, in D = 1, or 



-0.66, in D 



The energy per site of the two states including the 
quadratic terms are: £2-h = 2ejj — 4J 'D(en/ DJ) 2 , and 
e 2 _ V B - 2e VB - 4J(l/2)(-3/4) 2 where e V B = (—3/8) J 
is the energy per site of the VB state. In D = 1 the 
VB state wins over the Heisenberg state. In D = 2, 
the previous crude variational estimate would give the 
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Heisenberg state as favored. The VB state considered 
here is, however, very poor: its energy per site, neglect- 
ing the quartic term, is only — 0.375J, whereas it is well 
known that resonating VB states close in energy to the 
Heisenberg ground state can be constructed in D = 2.EII 

To show that the VB scenario is correct in ID, we cal- 
culate gaps and correlation functions using exact Lanc- 
zos diagonalizations of chains up to 14 sites, and Green 
Function Monte Carlo (GFMC) for longer chains employ 
ing the bias-correction scheme of Calandra and Sorella.Ej 
The trial wavefunction we have used is a Jastrow- 
type function (.Retrial) = {R\ exp{- J2(rr')[ v rr' (S? S*' + 
T r z T r z ,) - Aw rr ,{S z r S z r ,){T?T*)]}\R) where \R) is a con- 
figuration in the basis in which S z and T z are diago- 
nal. We allow, in general, long range tails in v rr i and 
w rr i , which are adjusted through a variational calcula- 
tion. Fig. [j] shows the finite-size value of the gap between 
the ground state and the first excited state both in the 
Sf ot = T z ot = 1 (□) and in the S z ot = T z ot = sectors 
(O)- The full symbols refer to -ffsT, whereas the open 
symbols refer to two decoupled Heisenberg chains (no 
quartic term in Eq. [l]) . All these excited states have mo- 
mentum 7r relative to the ground state. For the Heisen- 
berg case, the excited state is a triplet and the model 
has gapless excitations with gaps decreasing as 1/L for 
L — > oo (dashed line). For H^t, the situation is different: 
The lowest excited state is a singlet, with a finite size gap 
AE(S — 0) decreasing faster than 1/L for L — > oo. The 
triplet, instead, lies above with a gap AE(S — 1) extrap- 
olating to a finite value as L — > oo, AE(S = 1) w 0.5969. 
This is a clear signal of a degenerate infinite volume 
ground state with a spontaneous breaking of translational 
symmetry and a gap to all excitations, consistent with a 
VB phase. Further support to this interpretation can be 
gained by direct inspection of the spin-spin correlations, 
(SqS z ), and the dimer-dimer correlations along the chain 
((S§Sf)(S?S? +1 )). The inset of Fig. | shows a log-log 
size scaling of the peak of the two structure factors, lo- 
cated at momentum n. While the dimer-dimer structure 
factor peak, T ir (- 1 ) r (i S o s i)( S r S? +1 )) tends to diverge 
linearly with the length L of the chain, the usual spin- 
spin structure factor tends to a finite value as L — > oo. 

The possibility that a valence bond state might be the 
ground-state of the ST-model even in 2D was put forward 
in Ref.a. The argument, was based on a 5't/(n)-invariant 
model, due to Affleck,L3 which is unitarily equivalent to 
the ST-model in Eq. (|l|) for n = 4. It turns out that in 
the n — y oo limit such a S'{7(n)-invariant model should 
be characterized, in D = 2, by _a. plaquette resonating 
VB state,EJ as discussed in Ref.Efl. Such a phase has 
resonating VB plaquettes forming alternating disordered 
columns, but breaks translational invariance along the 
axis perpendicular to the disordered columns. 

We are now going to discuss the results we obtain with 
the GFMC for the ST-model on the square lattice. Fig. || 
shows the finite size gap between the ground state and the 
lowest excited triplet state for the ST-model on a L x L 
lattice as a function of the inverse volume. The data for 
the ST-model, at variance with those of the Heisenberg 



model, suggest the presence of a spin gap in the thermo- 
dynamic limit, although some care should be exercised 
because of the large error bar on the largest sample data, 
and the possibility of a non-trivial size scaling. The solid 
line is simply a straight line going through the two data 
points for a 6 x 6 and a 10 x 10 lattices. The extrapolated 
value of the spin gap is ~ 0.373 J. 

In principle, either a VB crystal with broken transla- 
tional symmetry or a homogeneous spin liquid is com- 
patible with the existence of a spin gap. In order 
to investigate more closely the possible kinds of long 
range order which might characterize the ground state 
of the ST-model in D = 2, we have calculated the 
expectation value of several spin-spin correlation func- 
tions by the forward walking technique.E3 Fig. || shows 
the results for the structure factors related to both 
standard spin-spin correlations ^2 r (— 1Y{SqS z ) (solid 
squares) and dimer-dimer correlations in x direction 
Er e lqr ((S§S z +£ ){S z S z +£ )) for different momenta: (0, tt) 
(circles), (ir, 0) (hexagons) (tt, tt) (triangles). For com- 
parison, we include the results for the (tt, tt) spin- 
spin structure factor of the D = 2 Heisenberg model, 
Er( _i ) r (^o^r) (open squares). We observe that all the 
dimer-dimer correlations considered show no sign what- 
soever of diverging linearly with the volume, i.e., we find 
that Sdim-dim(q) / L 2 — > 0. This finding rules out the pos- 
sibility of a crystalline VB state, even in the form of a 
plaquette resonating VB.EJ On the other hand, the stan- 
dard (tt, tt) spin-spin structure factor - signal of a Neel 
ordered state - , does not seem either to diverge linearly 
with the volume, i.e., S(tt, tt)/ L 2 also seems to go to zero. 
The latter finding is compatible with the presence of a 
gap suggested by the results in Fig. ^. Summarizing, al- 
though further work is necessary for a definite answer to 
the nature of the ST-model ground state on the square 
lattice, it seems fair to conclude that our data suggest a 
disordered VB state as the most likely candidate. 

In conclusion, we have studied a particular exchange 
Hamiltonian, Eq. (Q), which describes the low-energy 
physics of a Mott insulator with orbital degeneracy in 
the regime in which the inter-orbital singlet is the lowest- 
energy intermediate state available to virtual hopping. 
This regime tends to be further favored by a relatively 
strong electron-phonon coupling. The basic ingredients 
needed for a realistic system to be a candidate VB de- 
scribed in this work are: i) orbital degeneracy not triv- 
ially removed by cooperative JT and/or crystal field ef- 
fects; ii) relatively large molecules with a strong electron- 
phonon coupling, so as to make the inter-orbital singlet 
favored as compared to the Hund's rule triplet; iii) nar- 
row bands with relatively large on-site Hubbard U, so 
as to stabilize a Mott-Hubbard insulator at exactly one 
electron per site; iv) reduced dimensionality. At variance 
with the standard Hund's rule case (the high spin two- 
particle state has lowest energy), which would instead 
favor spin ferromagnetism, we find in our case a strong 
tendency to the formation of Valence Bond (VB) phases. 
In D = 1, the presence of a statically ordered crystal of 
VBs is firmly established. Given the presence of a size- 
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able spin-gap, weakly coupled chains should also have a 
crystal VB ground state. In D = 2 on the square lattice, 
a spin gap is suggested by our data, but no sign of crys- 
talline dimer order is found from the relevant correlation 
functions. A homogeneous liquid of resonating valence 
bonds appears as the natural candidate ground state for 
this model. 
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and ICTP was co-sponsored by INFM project HTSC, and 
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FIG. 2. Finite size gap between the ground state and the 
first excited triplet state for the ST-model (squares) and for 
the Heisenberg model (circles) on the square latticeEd The 
data for the ST-model are compatible with the existence of a 
spin gap in the thermodynamic limit. 

FIG. 3. Finite size structure factors obtained by GFMC 
with the forward walking technique. The important compo- 
nents of the Fourier transform of the dimer-dimer correlation 
functions look to grow more slowly than the volume. There- 
fore we can exclude the presence of any kind of valence bond 
solid. For the Neel order parameter, we find no clear sign of 
long range order. The size scaling of the Neel order parameter 
for the Heisenberg model is shown for comparison.^ 2 ! 
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FIGURE CAPTIONS 



FIG. 1. Finite size gaps for the lowest excited states of 
Hst in D = 1 (solid symbols) compared to the Heisenberg 
chain case (open squares). The excited state for the Heisen- 
berg case extrapolates to AE/J = iv 2 /(2L) + ■ • •, shown by 
the dashed line. The solid straight line, extrapolating to 
~ 0.5969 for L — > oo is obtained by a quadratic fit. In- 
set: Log-log plot of ( — l) r {SoSr) (open squares) and of 
£ r (-l) r ((SgSfHS^+i)) (solid squares). The dashed line 
has slope 1, for comparison. Top: Schematic illustration of 
the two crystal VB ground states of a linear chain. 
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